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We probe the dynamics of intermittent avalanches caused by steady addition of grains to a quasi- 
two dimensional heap. To characterize the time-dependent average avalanche flow speed v(t), we 
image the top free surface. To characterize the grain fluctuation speed Sv(t), we use Speckle- 
Visibility Spectroscopy. During an avalanche, we find that the fluctuation speed is approximately 
one-tenth the average flow speed, Sv » O.lw, and that these speeds are largest near the beginning 
of an event. We also find that the distribution of event durations is peaked, and that event sizes 
are correlated with the time interval since the end of the previous event. At high rates of grain 
addition, where successive avalanches merge into smooth continuous flow, the relationship between 



average and fluctuation speeds changes to 8v 



,1/2 



PACS numbers: 45.70.Ht, 83.70.Fn, 42.50.Ar, 45.70.Mg 
I. INTRODUCTION 

The behavior of granular media continues to pose 
significant challenges to both theory and experiment 
U 0, H, EL H| • One source of the richness in this field is 
that flows occur only after the forcing exceeds a thresh- 
old. The flow response is thus highly nonlinear, and can- 
not be understood as a small excitation above the static 
state. For example, grains in a hopper remain at rest 
unless the outlet underneath is greater than a few grain 
diameters. Grains on a vibrated plate remain at rest un- 
less the peak acceleration is greater than g. And grains 
on the surface of a pile remain at rest unless the slope is 
greater than an angle of maximum stability. In all cases, 
the resulting flows can be smooth and hydrodynamic- 
like at very high forcing, or temporally intermittent at 
low forcing just beyond threshold. 

In the case of surface flows, granular avalanches have 
been observed in a variety of geometries including a ro- 



tating drum d, 0, H, H, [n| j an inclined plane fill 
and a heap to which grains are added [1, M, MM 
Much effort concerns the size of the avalanches and the 
shape of the coarse-grained (a.k.a. hydrodynamic) ve- 
locity profiles. But what is the nature of the grain-scale 
dynamics? Fluctuations of individual grains away from 
the average velocity lead to inelastic collisions that dissi- 
pate energy and that ultimately control the nature and 
rate of the flows. While the microscopic grain dynamics 
is hence crucial to a fundamental understanding, it is par- 
ticularly difficult to measure in slow dense flows because 
the mean-free path and mean-free time for grain-grain 
collisions can be too short for high-speed video, which 
only captures surface behavior anyway. While diffusing- 
wave spectroscopy (DWS) is capable of measuring grain- 
scale fluctuations within the bulk, even for collision rates 
as high as 10 5 per second and for mean free paths as short 




it requires that the 
For time- 



as 10 -4 times the grain size 13, 
dynamics be time-independent 

dependent avalanching flows, DWS must be extended by 
consideration of higher-order intensity correlations be- 
cause the electric field statistics become non-Gaussian 
fl5l fl6| . However, interpretation of such data relies on 
the assumption that flows start and stop instantaneously, 
such that the average and fluctuation speeds remain con- 
stant during the course of an event. 

In this paper we apply a time-resolved multi-speckle 
dynamic light scattering technique known as speckle- 
visibility spectroscopy (SVS) [22l |23|| to avalanches on a 
granular heap confined between two parallel plates. This 
method is applicable to slow dense flows, just like DWS 
and the higher-order extensions, but is capable of resolv- 
ing the evolution of dynamics throughout the course of 
an avalanche event. Thus we measure in detail not just 
the statistics of avalanche sizes, but also their dynamics. 
And by comparison of the resulting fluctuation speeds 
with coarse-grained average speeds found by high-speed 
video, we find that the nature of microscopic dynamics 
is different during an avalanche than during continuous 
flow. 
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II. EXPERIMENTAL METHODS 

A. Granular system 

A quasi-two-dimensional heap is created by steady ad- 
dition of grains between two parallel vertical walls, made 
of static-dissipating Lucite plates closed at bottom and 
along one vertical edge, as in Refs. 0,13. The gran- 
ular medium is dry, cohesionless glass beads with diam- 
eter range 0.25 — 0.35 mm, repose angle 8 r = 25°, and 
density p = 1.5 g/cm 3 . The inner wall dimensions are 
28.5 x 28.5 cm 2 , so that the length along the slope of the 
heap is 28.5 cm/ cos9 r = 31 cm. The channel width is 
w = 9.5 mm, equal to approximately thirty grains across. 
The flux Q of grains onto the heap may be varied widely 
in increments of 0.005 g/s, and be held constant for ex- 
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tended periods, as follows. A large funnel is positioned 
above the back of the channel, filled with grains, and al- 
lowed to drain under gravity. The output is connected 
to a valve, which uses a moveable knife-edge to divert 
the desired flux onto the heap and to expunge the rest 
into a storage bucket for manual recycling. The free- 
fall of grains onto the top of the heap is broken by an 
Aluminum bar, sandwiched between the plates to form a 
funnel with 1.5 cm outlet 3 cm above the top of the heap. 



Several distinct flow regimes are observed depending 
on the flux Q of grains. Above a critical flow rate of 
0-36 g/s, 



Q c = 0.36 g/s, the same value as in Ref. [15j, the flow 
along the surface of the heap is smooth and continuous 
[HI, US HH . Below Qc, the flow is intermittent, accom- 
plished by a series of discrete avalanches that start at the 
top of the heap. About every third or fourth avalanche 
is a large, system-wide event in which the entire surface 
flows and expunges sand out the bottom of the channel; 
after such an event, the surface of the heap is smooth 
with a constant slope along its entire length. Other 
avalanches are smaller and stop before reaching the bot- 
tom; after these events, the surface of the heap is un- 
even, with discrete steps between long regions of constant 
slope. Very far below Q Cl the behavior is quasi-static in 
that successive avalanches are independent of one another 
and also in that variation of Q affects the time between 
events but not the dynamics during flow. As Q c is ap- 
proached from below, successive avalanches merge to an 
extent set by the value Q c — Q. For the avalanche exper- 
iment reported below, the flux of grains is held fixed at 
Q = 0.07 g/s. This lies in the quasi-static regime, but 
only slightly below the point at which successive events 
start to merge in order to maximize the number of events 
observed per unit runtime. At the very end of the pa- 
per, grain dynamics during avalanche are compared with 
those under continuous flow at grain fluxes across the 
range Q c — 0.36 g/s < Q < 3.3 g/s. 



Before describing the diagnostic tools and their ap- 
plication to behavior near the top free surface, we first 
remark upon the three-dimensional character of the ve- 
locity profiles in the continuous flow regime. Using parti- 
cle imaging velocimetry, as described below, we find that 
the flow speed at the side walls decreases nearly expo- 
nentially with depth. The decay length is comparable to 
the channel width. We also find that the flow profile is 
approximately parabolic across the top free surface, but 
with wall slip such that the speed at the center is 1.4 
times faster than the speed at the wall. For a known 
flux of 2.5 g/s, multiplying these forms and integrating 
over the cross section of the heap gives a flux estimate 
of 2.2 g/s. The good agreement suggests that flow along 
the top and side surfaces is representative of flow within 
the bulk. 



B. Diagnostics 

Next we describe technical details for two methods for 
characterizing time- resolved avalanche dynamics. Both 
are applied to the top surface of the granular heap, over 
the region between 2 and 4 cm from the outlet as mea- 
sured along the slope. The only avalanches that pass 
through this field of view are large ones, which span 
the entire surface of the heap. The first diagnostic is 
a variant of Particle-Image Velocimetry (PIV) [H |28| . 
which gives v(t) the time-dependent coarse-grained hy- 
drodynamic surface flow speed. While originally devel- 
oped for fluids, this method is well suited for granular 
media, as in eg. Refs. [HIM IM HI- The second diag- 
nostic is Speckle- Visibility Spectroscopy (SVS) [22l |23|. 
which gives Sv(t) the time-dependent speed of micro- 
scopic grain-scale fluctuations around the hydrodynamic 
flow. These methods are implemented one at a time us- 
ing an 8 bit Basler digital linescan camera, with 2 x 1024 
pixels and 58 kHz maximum frame rate, controlled us- 
ing Lab VIEW v7.1 and the National Instruments Vi- 
sion Toolkit. Images are continuously captured and pro- 
cessed, so that raw video data need not be saved and so 
that extremely long duration runs are possible. 



1. Particle imaging 

The time-dependent average speed, v(t), of grains at 
the top free surface may be deduced from the spatial 
cross-correlation of successive images, as follows. Here, 
a bright halogen lamp is placed about 1 m away from 
the heap, shining down between the Lucite plates. The 
linescan camera is placed about 20 cm away, with optical 
axis normal to the heap and and with the two rows of 
pixels oriented parallel to the flow direction. The camera 
is fitted with a lens, such that the field of view is a 2 cm 
long strip, 39 ^m wide, located half way between the side 
walls. Under these conditions, surface beads reflect light 
back to the camera and appear as a dark circles with 
central bright spots with object size of about 1/4 bead 
diameter and image size of size 4 pixels. The frame rate 
of the camera is set to 1 kHz, such that for a typical 
flow speed of 4 cm/s the beads move about 2.5 pixels per 
frame. 

While beads are thus imaged, it is not necessary to 
identify and track their positions individually. Rather, 
the ensemble average speed of all imaged beads is found 
by the average displacement between successive frames. 
This displacement is computed as the peak position of 
the spatial cross-correlation of successive images. Here, 
correlations are found by Fourier methods, and the peak 
positions are identified by parabola fit to the cross- 
correlation function. 

Example image data and flow speeds for a typical 
avalanche are shown in Fig. [TJ including a blow-up of 
the beginning and end of the event. The top row shows 
space-time plots of raw grayscale images. Before and 
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FIG. 1: (Color online) Time-stacked linescan images of the surface of the heap for an avalanche (a) and for zoom-ups of it 
turning on (b) and turning off (c). Sand grains appear as bright lines with slope proportional to their speed down the channel. 
The slope is measured by computing the maximum of the cross-correlation function for consecutive time snapshots. The bottom 
row shows the flow speeds extracted from the corresponding images above. 



after the event, when grain are at rest, the individual 
bright spots all remain at the same location and hence 
cause horizontal streaks in the plot. During the event, 
the bright spots move and hence cause streaks with slope 
set by individual grain speeds. The average streak slope, 
computed as per above, gives the ensemble-average grain 
speed shown in the bottom row. Since grains diffuse lat- 
erally during flow, individual streaks last for only a finite 
time duration. 

Several points are to be noted in Fig.[TJ First, the grain 
speeds are not constant during the avalanche. Rather, 
v(t) increases to a maximum at the beginning of the event 
and then gradually decreases back to zero. A front of 
flowing sand is observed to sweep through the field of 
view in part (b); this sets the rise time of v(t) rather 
than the acceleration of grains from rest. By contrast, 
grains in the field of view in part (c) all appear to come 
to rest at the same time. The final approach to rest is 
very abrupt, nearly but not quite discontinuous. The 
scale of noise in v(t) data is consistent with the estimate 
Av = f/{MyN) — 0.5 cm/s, where / is the frame rate, 
M is the magnification in pixels per cm, and N is the 
number of beads in the field of view. 



2. Speckle visibility 

The time-dependent fluctuation speed, Sv(t), of grains 
at the top free surface may be deduced from the the 
visibility of laser speckle formed by backscattered light. 
This method has been dubbed Speckle- Visibility Spec- 
troscopy (SVS) [23j], and has been applied to grains sub- 
ject to periodic vibration [22j], to colloids after cessation 
of shear [33[, and to foams subject to coarsening [23| . 
Closely related methods are laser-speckle photography 



[3J and Time-Resolved Correlation (TRC) [3J, [3J, [3J . 
Here, coherent light from a Nd:YAG laser, wavelength 
A = 532 nm, power 4 W, is expanded to a Gaussian di- 
ameter of 1.3 cm and is directed normally onto the heap 
at the same location that average speeds were observed. 
An aperture blocks the beam near the channel walls, so 
that the illumination spot is roughly 0.95 x 1.3 cm 2 . In- 
cident photons diffuse in the medium with a transport 
mean free path of a few grains. Therefore the typical 
photon emerges from the sample after a couple scattering 
events from grains close to the surface [38| . The linescan 
camera is placed about 15 cm away, with optical axis nor- 
mal to the heap. Now, however, the lens is removed and 
is replaced by a laser line filter to eliminate room light. 
Under these non-imaging conditions, the backscattered 
produces a speckle pattern in the plane of the camera 
such that the speckle size is about 0.7 pixels. As grains 
move relative to one another, the speckle pattern fluctu- 
ates and hence appears visible only if the camera frame 
rate is fast compared to the speckle lifetime, which is set 
roughly by the time X/Sv for adjacent scattering sites to 
move one wavelength apart. Note that for uniform trans- 
lation of the sample, with no motion of grains relative 
to one another, the speckle pattern changes much more 
slowly over the time L/v required for a new ensemble of 
grains to come into the field of view L. Hence the idea 
of SVS is to deduce the grain fluctuation speed Sv from 
the visibility of the speckle for a given exposure duration 
T. Here the camera is operated at maximum frame rate, 
58 kHz, giving T — 17.24 /is. For convenience, only the 
central 2 x 430 pixels are processed. The laser intensity 
is adjusted so that the average grayscale level is 50. 

The visibility of the speckle may be quantified by the 
variance of intensity levels, V(T) cx (I 2 )t — (I) 2 , where 
(• • -}t denotes the average over pixels exposed for dura- 
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FIG. 2: (Color online) Time-stacked line scan images of the speckle pattern for an avalanche (a) and zoom-ups of it turning on 
(b) and turning off (c). The middle row shows the corresponding variance ratios, which are transformed into the fluctuation 
speeds in the bottom row. 



tion T. Note that the average intensity is independent 
of T, so no subscript is placed on (I). The proportional- 
ity constant of V(T) is set by the laser intensity and the 
ratio of speckle to pixel size. It may be neatly canceled 
by considering the variance ratio V(2T)/V(T), where the 
numerator is found from a "synthetic exposure" equal to 
the sum of successive images. For diffusely backscattered 
light from particles moving with random ballistic motion, 
the theory of SVS [22|, [23| gives the variance ratio as 



V(2T) _ er 
V(T) ~ T[er 



- 1 +Ax 



l + 2x) 



(1) 



where x = (Air5v/\)T. This equation may be numer- 
ically solved for Sv vs time in terms of data for vari- 
ance vs time. A rational approximation that has the 
correct small and large x limits of 1 — 2x/3 + 0(x 2 ) and 
(4 + l/x)/8 + 0(1/ x 2 ), respectively, and that may be in- 
verted analytically for a good initial guess, is given by 
V{2T)/V{T) w (1 + 7x/& + x 2 )/(l + llx/6 + 2x 2 ). 

Example speckle images and SVS analysis for a typi- 
cal avalanche are shown in Fig. [2] including a blow-up of 
the beginning and end of the event. The top row shows 
space-time plots of raw grayscale images of the speckle 
pattern. Before and after the event, when grain are at 
rest, the speckles are nearly static and hence appear as 
horizontal streaks. The variance of grayscale levels is 
nearly maximum, both as measured over exposure times 
T and 2T . Hence the variance ratio is nearly one and 



the random grains speeds are nearly zero, as shown in 
the second and third rows. However the speckle streaks 
do not extend indefinitely in the top row of Fig. [21 with 
speckles lasting on the order of a few seconds. Over this 
time the grains experience wavelength-scale motion, due 
perhaps to ambient vibration, to the addition of grains 
at the top of the heap, or even to thermal expansion 
and contraction [33 |. The motion during the avalanche 
is much faster, so that the random bright-dark pattern of 
speckles becomes a washed-out blur readily visible only 
near the beginning and end of the event. Therefore the 
variance of grayscale levels decreases, more so for increas- 
ing exposure duration, as quantified by the variance ratio 
V(2T) /V(T) displayed in the middle row of Fig. EH Using 
Eq. Jl} , this data gives the random fluctuation speeds Sv 
displayed in the bottom row. Note that, like the average 
speed v, the fluctuation speed is maximum at the begin- 
ning of the event, then it decreases gradually and finally 
halts abruptly. Also note that the two speed scales differ 
by a factor of roughly ten, Sv = O(0. l)v. 



III. RESULTS 

Example ten-minute time-traces of the average flow 
speed v(t), from imaging, and of the fluctuation speed 
Sv(t), from SVS, are displayed in Fig.[3J Each spike rep- 
resents one avalanche, in which both speeds suddenly rise 
from zero, persist briefly, and then vanish. Note that sue- 
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FIG. 3: (Color online) (a) Average and (b) fluctuation speed 
timetraces for independent observations with the same sand 
flux Q = 0.07 g/s. The on and off times, T on and T a s, for the 
n th avalanche event are defined as shown. 



cessive avalanches are roughly similar in size, duration, 
and spacing. In this section we analyze two 37 hour runs, 
spanning 1436 events in v(t) and 1472 events in Sv(t). We 
find that the root-mean-squared speeds during flow are 
sj (v 2 ) = 3.0 cm/s and -J (Sv 2 ) = 0.37 cm/s, and that the 
fraction of time spent in flow and at rest are /i = 0.2 and 
/o = 0.8, respectively. First we consider the times T on 
and T ff, defined in Fig. [3l for the duration of events and 
for the intervening quiescent periods. Then we consider 
switching functions that describe the probability for the 
heap to change between flow and rest states after a given 
time delay. Lastly we consider the dynamics of an av- 
erage event and the relation between the instantaneous 
average and fluctuation speeds. 



A. Event durations 

Histograms for the rest time T a g between avalanches, 
and for the duration T on of events, are collected in 
Figs. 2h.-b- Whether measured by imaging or SVS, the 
results show that the distribution of rest times is strongly 
peaked with an average of (T Q ff) = 61 s and a standard 
deviation of 34 s. The minimum rest time is about 10 s, 
which is 1.5 standard deviations below the average; there- 
fore successive avalanches are well-separated and quasi- 
static. As for the rest times, the distribution of event 
durations is also strongly peaked, with an average of 
(Ton) = 16 s and a standard deviation of 4 s. While 
both distributions are skewed toward longer times, they 
have a rapid final decay inconsistent with a power-law. 
Hence the events are quasi-periodic; the average period 
is <T off + T on ) = 77 s. 

A length scale quantifying event size may be defined 
by integrating the flow speed over the event duration, 
L = J v(t)dt. The distribution of event lengths, shown 
in Fig. |4j:, is once again peaked with an average of 
(L) = 45 cm and a standard deviation of 14 cm. The his- 
tograms for L and T on have slightly different shapes be- 
cause v(t) decreases throughout the event. Note that (L) 



is about one standard deviation longer than the distance 
along the slope of the heap, consistent with avalanches 
that sweep through the entire system. Further note that 
the average thickness of the flowing layer is given by 
mass conservation as Q(T on +T a s) / (p{L}w), which equals 
about three grain diameters - consistent with visual ob- 
servation and thinning of the layer during the course of 
an event. 

Next we consider possible correlations between the rest 
times and event sizes. As in the example time trace, 
Fig. [3l the event index n is defined such that T g(n) 
is the rest time immediately prior to the n th avalanche, 
which has duration T on (n). Scatter plots of flow and 
rest times vs preceding rest time, and of rest and flow 
times vs preceding flow time, are displayed in Fig. [5^-d. 
The only noticeable correlation is in part-c, which shows 
that T on (n) grows with increasing T g (n). Intuitively, for 
longer rest times more grains accumulate at the top of the 
heap and so the next event is larger. There appears to be 
no such correlation between successive events. Thus T Q g 
is determined at random from the distribution shown in 
Fig.H but T on is not. 



B. Switching probabilities 

To study the dynamics by which avalanches start and 
stop, we begin by computing the autocorrelation of the 
v(t) and Sv(t) time trace data. The results, plotted in 
Figs. [6^,-b, display a linear decay followed by damped 
oscillations to a constant. This is consistent with quasi- 
periodic behavior. Note that the initial decay time is 
set by (T on ) and that the first peak occurs at the pe- 
riod (T ff + T on ). The detailed shape of the speed au- 
tocorrelations is determined both by the variation of 
the speed during an event, as well as by the statistics 
by which the flow starts and stops. To investigate the 
relative importance of these contributions, we compare 
with expectation based on telegraph-approximated sig- 
nals v(t) w \J (v 2 )x{t) and Sv(t) w \/ (5v 2 )x(t), where 
x(t) is equal to during rest and is equal to 1 dur- 
ing flow. The resulting autocorrelations, displayed in 
Figs.[6K-b, have nearly the same shape as the actual speed 
autocorrelations. Note that the initial and final expec- 
tations are set by the mean-squared speeds multiplied 
by (x 2 ) = fi and (x) 2 — fi 2 , respectively. The good 
agreement between actual and telegraph-approximated 
autocorrelations implies that the variation of speed dur- 
ing an event has only minor consequence. This supports 
the validity of the analysis of the higher-order intensity 
correlation data presented in Refs. [la, EH]. 

To fully characterize the dynamics of switching be- 
tween flow and rest states, we employ standard prob- 
ability functions as in Refs. [HI, [l6[. Thus we define 
Pij (t) as the conditional probability for the system to be 
in state j at time t + r if it started in state i at time 
t. By convention subscript denotes a state of rest, and 
subscript 1 denotes at state of flow. All four of these 
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FIG. 4: (Color online) Distributions of (a) off times, (b) on times, and (c) lengths L — J t ' o// v(t)dt. Blue curves correspond to 
results based on v(t) from particle imaging experiments, in which 1436 events were observed. Red curves correspond to Sv(t) 
from speckle-visibility experiments, in which 1472 events were observed. The averages are {T g)g v = 60.0 s, {T g) v = 62.3 s, 
(T on ) Sv = 14.9 s, (T on ) v = 16.6 s, and (L) = 44.8 cm. 




FIG. 5: (Color online) Scatter plots of on and off times for successive events, as labeled. As defined in Fig. the temporal 
sequence is T ff(n), T on (n), T s(n + 1), T on (n + 1), etc. Significant correlation exists only in part-c, as indicated by the solid 
gray line. 



inter-related functions may be computed from the speed 
vs time data, in terms of the autocorrelation (x(O)x(r)) of 
the corresponding telegraph signals. The relevant iden- 
tities are (x(O)x(t)) = /iP u (t), P 00 (r) + P 01 (t) = 1, 
Pio{t) + Ph(t) = 1, and / P i(t) = /iPio(r). In addi- 
tion we define two more functions, Pq(t) and Pl(t), as 
the probabilities to be in the same on or off state at time 
t + t as at time t, with no changes of state in between. 
These may be computed by averaging over all off times as 
P (r) = ('[1 - T/T oB (n)]H(T oS (n) - r)), and similarly for 
Pi(t), where H(t) is the Heaviside function. The func- 
tion Po(t) is particularly crucial for analysis of dynamic 
light scattering data O EH . 

Results for the switching probabilities are collected in 
Fig. [7] The initial decays of Po (r) and Poo (r) are both 
linear, 1 — r/T g. While the former decays fully to zero, 



since all rest states have finite duration, the latter os- 
cillates and decays to /o = 0.8 due to contributions 
from other rest states. Note that the avalanches are 
quasiperiodic, in that about three oscillations occur be- 
fore the full decay to /o- Similar statements hold for the 
analogous flowing state switching functions. The func- 
tional form for Pq (r) is similar to that found in Ref. [l6| 
using multiple light scattering. Namely, it is faster 
than exponential and well-described by [1 + 3(r/T g) + 
4(r/T off ) 2 + (8/3) (r/T off ) 3 ] exp(-4r/T off ). This function 
is plotted as a solid gray curve; its limiting behavior is 
1 - (r/T off ) + (32/15)(r/T off ) 5 + O(r/T off ) 6 . 
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FIG. 6: (Color online) Autocorrelation functions of (a) 
v(t) and 5v(t) data. The average values (Ton), (Toff), and 
/i = ET on /(£T on + ET ff) are computed from distributions 
in Fig. [4] Solid curves are autocorrelations of actual data, 
while dashed curves are autocorrelations of the telegraph- 
approximated signals. 
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FIG. 7: (Color online) Probability functions for switching be- 
tween states vs delay time, as labeled. Blue and red curves 
correspond to results based on v(t) and Sv(t) data, respec- 
tively. The solid gray curve is [1 + 3(r/T off ) + 4(r/T off ) 2 + 
(8/3)(r/T off ) 3 ] exp(-4r/T off ). 



C. Average event dynamics 



for v(t) and 5v(t) decrease nearly linearly with time dur- 
ing a span of about ten seconds. This is indicated by the 
dashed black curves, which trace along the crest of the 
probability distributions. After 10-15 s, the typical be- 
havior becomes less well-defined since some events stop 
and some continue. In parts (b,d), the typical behavior 
at the end of an event is revealed by binning the speed 
data for all events vs the time t — t Q g until the avalanche 
ceases. As before, the crest of the distributions is indi- 
cated by a black curve. The data show that the speeds 
continue the initial linear decrease with time up until 
about 5 s before cessation. During the final stage, both 
v(t) and 5v(t) decrease ever more rapidly with time until 
reaching zero. While there is no discontinuity in speeds 
at t — tag, the slopes vanish abruptly. 

To examine the relationship between flow and fluc- 
tuation speeds, we display a log-log parametric plot of 
(Sv(t)) vs (v(t)) in Fig. H' Since PIV and SVS data were 
collected separately with the same camera, we do not 
have simultaneous flow and fluctuation speed data for 
the same avalanches. Therefore, the results on display 
are for a typical event given by the black curves in Fig. [8] 
tracing along the crests of the probability distributions. 
In other words, (. . .) represents the average taken over an 
ensemble of events. When events are aligned according to 
start time, t—t on , the relation between (Sv(t)) vs (v(t)) is 
displayed as a orange dashed curve. And when events are 
aligned according to stop time, t — t ff, the relation be- 
tween (Sv(t)} vs (v(t)} is displayed as a solid purple curve. 
These two curves agree well, but the latter has a larger 
range since both speeds vanish at the end of the event. 
The relation is approximately (Sv) — 0.1(w), as shown 
by the dashed line with slope 1; during an an avalanche 
event, the fluctuation speed is proportional to the flow 
speed, but is ten times smaller. For smooth continuous 
flow at high flux, Q > Q c , a different relation is found. 
The data, shown by open green circles in Fig. [9l are con- 
sistent with (Sv) ~ \J (v) as shown by the dashed line 
with slope 1/2. Similar behavior was observed in hopper 
flow, where the exponent was approximately 2 /3 [171 ] . For 
continuous flows, the fluctuation speed decreased more 
gradually than the average flow speed as the flux is low- 
ered. This results in proportionately greater dissipation, 
and hence a transition to intermittent flow at nonzero 

Qc- 



In this final section we consider the dynamics of indi- 
vidual avalanches, and how the average and fluctuation 
speeds evolve during the course of an event. Since there 
is a range of avalanche sizes and durations, the results 
are shown in Fig. [8] as probability densities for there to 
a given speed at a given time. In parts (a,c), the typical 
behavior at the beginning of an event is revealed by bin- 
ning the speed data for all events vs the time t — t on since 
the start of flow. The initial rise from zero is quick, set 
by the sweep of a moving front through the field of view. 
Once the speeds reach a maximum, the average results 



It is curious to note in Fig. [5] that (8v) vs (v) data for 
continuous and intermittent regimes are not disjoint. At 
the beginning of an avalanche, and for continuous flow at 
fluxes slightly greater than Q c , the average flow speeds 
can be the same. Across this overlap, the fluctuation 
speeds in the two regimes are in good agreement. There- 
fore, the nature of granular heap flow at the beginning 
of an avalanche is remarkably similar to that for contin- 
uous flow at low flux. In this sense, there is a smooth 
crossover from (Sv) ~ (v) for avalanches to (Sv) ~ v/ (v) 
for continuous flow. 
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FIG. 8: (Color online) Probability density maps (rainbow colors), and most probable avalanche dynamics (black curves), for 
the average speed v(t) and fluctuation speed 5v(t) vs time t during an avalanche. All observed avalanches are included by lining 
up their individual time traces according to either when they turn on (a,c) or else according to when they turn off (b,d). The 
bin sizes are 0.02 cm/s x 0.1 s for v(t) data and 0.002 cm/s x 1.72 ms for Sv(t) data. 
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FIG. 9: (Color online) The relation between (Sv) and (v) for 
both continuous flows (open green circles) and avalanching 
flows (curves). The dashed orange curve represents {Sv(t — 
ton)) vs (v(t — ton)), based on the dashed curves in Fig. [HJl 
and c. The solid purple curve represents {Sv(t — t ff)) vs 
(v(t — toff)), based on the solid curves in Fig. [8}j and d. 



IV. CONCLUSION 

Speckle- Visibility Spectroscopy has permitted us to ob- 
serve the instantaneous velocity fluctuations of sand par- 
ticle for continuous and avalanching flows over a large 
dynamical range. Particle-image velocimetry has allowed 
us to observe the instantaneous flow velocity of sand par- 



ticles in continuous and avalanching flows over a simi- 
larly large dynamical range. Together, these observation 
methods provide Sv and v measurements that allow us to 
significantly improve upon previous studies [TEl . [Ifjj and 
to thoroughly characterize both the macroscopic and mi- 
croscopic dynamics of avalanches. 

The ability to observe thousands of distinct avalanches 
over the course of 74 hours of observation has permitted 
us to uncover the full time-on and time-off distributions. 
We learn that even for constant sand addition flow rates, 
avalanches come in a wide variety of sizes, and the on 
and off times of avalanches are sampled from distribu- 
tions that are peaked but non-Gaussian. The avalanches 
are quasi-periodic over only a few avalanche cycles, and 
the flow periodicity is essentially decorrelated after three 
cycles. 

The ability to fully observe the instantaneous flow and 
fluctuation dynamics of independent avalanches has al- 
lowed us to uncover their rich and characteristic dynam- 
ical shape. There is a sharp wavefront at the head of 
the avalanche, continuous deceleration in the middle, and 
abrupt cessation of flow at the end. Thus, even though 
the simple square wave of the telegraph model works 
well on long timescales, it is a poor approximation to 
the richly detailed shape of individual avalanches in the 
flowing state. 

Finally, because we measure Sv and v instantaneously 
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for tens of hours and thousands of avalanches, we are able 
to directly test the validity and range of applicability 
of previously reported functional forms. The result is 
that we learn that the functional form Sv ~ v 1 I 2 is only 
correct for highly fluidizcd continuous granular flows. For 
intermittent flows, this law progressively fails and gives 
way to Sv ~ v precisely at the critical flow rate crossover 
Q c that separates continuous from avalanching flows. 

These results combine to give a detailed picture of 
avalanches at both the microscopic and macroscopic lev- 
els and with the range necessary to fully characterize 
dynamics. The microscopic dynamics of sand particles 
depend on the macroscopic flow state of the pile. For 
example, the "granular temperature" in a continuously 
decelerating avalanche is lower than would be expected 
by extrapolating back from the continuous flow relation 
Sv ~ v 1 / 2 , in which correlated particle collisions are a 
necessary ingredient [l7j . Instead, microscopic dynam- 



ics of intermittent avalanches fit far better to the form 
Sv ~ v. New theoretical insights may uncover how chang- 
ing flow-structure of decelerating avalanches [l(| affects 
the collision rules of sand particles at the microscopic 
level, which are empirically found to be different on ei- 
ther side of the critical transition point Q c - 
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